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© Method and apparatus for automatically producing representations of three-dimensional horizons from processed 
seismic data. 

© An electronic system automatically produces repre- 
sentations of three-dimensional recommended horizons 
(50a-50e) from processed seismic data (40), Significant 
events, in the form of turnings, are identified in all of the digit- 
ized traces that make up the processed seismic data. Pairs of 
significant events in adjacent traces are determined. The 
system determines these pairs through a mutual nearest 
neighbor technique, together with a technique for detecting 
similar signs. The system measures the curvature of the path 
formed by the paired turnings in successive digital traces and 
terminates the path at any point where the curvature exceeds 
a predetermined limit. The system then looks for three- 
dimensional continuity by looking for a closed loop of at least 
four turnings making up four pairs of turnings. These tech* 
niques yield at least one horizon represented in a first format. 
The recommended horizon <50a-50e), represented in a second 
format, is essentially extracted by the system through a recur- 
sive 'walking" procedure for identifying all of the turnings In 
the recommended horizon and designating those turnings in 
the second format. At the same time, vertical uniqueness of 
each of the turnings inthe recommended horizon is examined 
and if a subsequent turning is found to be on the same digital 
trace as a prior detected turning, the subsequent turning is not 
permitted. The second format is used for storage of the hor- 
izon and subsequent display. 
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seismic data. 

© An electronic system automatically produces represen- 
tations of three-dimensional recommended horizons (50a- 
50e) from processed seismic data (40), Significant events, in 
the form of turnings, are identified in all of the digitized 
traces that make up the processed seismic data. Pairs of 
significant events in adjacent traces are determined. The 
system determines these pairs through a mutual nearest 
neighbor technique, together with a technique for detecting 
similar signs. The system measures the curvature of the path 
formed by the paired turnings in successive digital traces 
and terminates the path at any point where the curvature 
exceeds a predetermined limit. The system then looks for 
three-dimensional continuity by looking for a closed loop of 
at least four turnings making up four pairs of turnings. These 
techniques yield at least one horizon represented in a first 
format. The recommended horizon (50a-50e), represented in 
a second format, is essentially extracted by the system 
through a recursive "walking" procedure for identifying all of 
the turnings in the recommended horizon and designating 
those turnings in the second format. At the same time, 
vertical uniqueness of each of the turnings in the recom- 
mended horizon is examined and if a subsequent turning is 
found to be on the same digital trace as a prior detected 



turning, the subsequent turning is not permitted. The second 
format is used for storage of the horizon end subsequent 
display. 
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METHOD ANO APPARATUS FOR AUTOMATICALLY PRODUCING REPRESENTATIONS 
OF THREE-DIMENSIONAL HORIZONS FROM PROCESSED SEISMIC DATA 

RAPKfiRfHJND OF THE INVENTION 
Field of the Invention 

This invention relates to the analysis of seismic data and 
specifically to an electronic system for automatically producing 
5 horizons from processed seismic data. 

DESCRIPTION OF THE PRIOR ART 

In the past, human interpreters solved the problem of 
producing horizons of three-dimensional seismic events mentally, 
using paper and colored pencils, or computer graphics, to draw a 
10 seismic event as it intersects successive two-dimensional seismic 
sections of a data volume. Once two-dimensional events along the 
three-dimensional surface have been picked, the two-way times to 
each of the 2-D events are determined and a variety of different 
types of maps are drawn to illustrate information about the event 
15 surface. The primary disadvantage of this approach is that it is 
very time-consuming and tedious for the human interpreters. The 
task of drawing maps or horizons of event surfaces is easily 
automated once the events are picked, but the picking is very 
time-consuming. A prior art method for picking and mapping 
20 events in 2-D seismic sections has been achieved. This previous 
approach to determining significant changes by amplitude 
zero-crossings in seismic traces essentially amounts to an 
assertion that the significant changes of interest are primarily 
determined by points where the trace changes sign. This technique 
25 is not well-suited to mapping seismic events whose amplitudes 
change signs. 



BRIEF SUMMARY OF THE INVENTION 

A digital computer system 
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representations of three-dimensional recommended horizons from 
processed seismic data. The seismic data is in the -form of a 
3-D volume of individual seismic traces in a digital format. 

Significant seismic events, in the form of turnings, are 
identified in all of the digitized traces that make up the 
processed seismic data. These turnings are defined as 
zero-crossings in the Nth derivative of a seismic trace. The 
digitized seismic traces of this Invention are represented by a 
series of points. An approximation of the Nth derivative at each 
point is made. A zero-crossing in the Nth derivative determines 
a turning point. A turning occurs at the location of each such 
zero-crossing and the sign of the turning is the sign of the 
(N+l)th derivative at that point. This technique is described 
and claimed in copending U. S. Patent Application, Serial No. 
(Tl-1 0060 1 4tq,iTo*f-. 

Next, in the X and Y directions, through the 3-D data 
volume, pairs of turnings are determined. In this preferred 
embodiment, pairs of turnings are determined through the use of 
the mutual nearest neighbor technique. In this technique, a pair 
occurs wherever two turnings, taking one each from adjacent 
traces, are mutual nearest neighbors. That is, point A is the 
nearest neighbor of point B, and point B is the nearest neighbor 
of point A. This mutuality presumes that the two turnings 
occurred on the same three-dimensional seismic surface (horizon). 

The signs of the turnings are then determined, and if alike, 
further assure that points A and B are in the same horizon. The 
sign of each turning is the sign of the (N+l)th derivative at 
each zero-crossing. The determination of pairs of turnings is 
described and claimed in copending U. 5. Patent Application, 
Serial No. ( TI- 10058) bfcq , GoZ. 

The curvature of the path formed by the paired turnings in 
successive digital traces, in both the X and Y direction, is 
examined. If the curvature is too extreme, the path is 
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terminated. Restricting pairs of turnings in this manner 
eliminates some implausible 3-D seismic events. 

Three-dimensional continuity is then tested. This is done 
by looking for closed loops containing four turnings formed by 

5 four pairs of turnings. The existence of such a closed loop 
helps to assure that the seismic event has three-dimensional 
extent as a surface, and -not just as a curve. When a pair is not 
a part of such a closed loop, the pair is eliminated. This 
examination and modification of tentative horizons is described 

10 and claimed in copending U.S.. Patent A PP 1 i cat i on Serial No. 

( TI-10064) bt><»£03. 

By reason of the system performing the above steps, one or 
more horizons have been defined within the three-dimensional 
block of data. 

15 The defined horizons may now be, in essence, extracted from 

the 3D data volume for ease of display and handling in general. 
This is done by the system through reformatting from the network 
format, resulting from the above treatment, to an array format. 
At the same time, another analysis is performed to determine 
20 that defined horizons are not interconnected. A test for 
vertical uniqueness is performed by excluding any subsequent 
turning on the same digital trace as a prior detected turning. 
This vertical uniqueness assures that a horizon cannot spiral. 
An array format provides a more convenient representation of each 
25 horizon. The arrays are better suited for storage and display 
than the networks. This extraction of horizons is described and 
claimed in copending U. S. Patent Application. Serial No. 

(Tl-10059) • 

The object of this invention is to provide a system that 
30 automatically defines three-dimensional horizons from a 
three-dimensional volume of processed seismic data. This and 
other objects wi 1 1 be made evident in the detailed description 
that fol lows. 
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RR1EF DESCRIPTION OF THE DRAWINGS 

Figure 1A illustrates a three-dimensional volume of seismic 

data . 

Figure IB illustrates three-dimensional horizons produced 
from the volume of Figure I A, using the method and apparatus of 
this invention. 

Figure 2 is a block diagram of the computer system of this 
i n vent ion. 

Figure 3 illustrates a typical trace is analog form. 

Figure 4 illustrates two analog traces in the application 
of the mutual nearest neighbor criterion. 

Figure 5 is a reproduction of Figure 4 except that the 
effect of also requiring like signs is illustrated. 

Figures 6A-6H illustrate traces and indicate horizons 
determined from zeros of successive derivatives taken with N=0 at 
Fig. 6A and N=7 at Fig. 6H. 

Figure 7 illustrates the low curvature criterion. 

Figure 8A illustrates 3-D continuity of four turnings 
forming four pairs of turnings. 

Figure 8B, from a top view, illustrates 3-D continuity. 

Figure 9 illustrates vertical uniqueness. 

Figure 10 illustrates transitivity of the pair relation to 
reformat to the array format. 

DETAILED DESCRIPTION OF THE INVENTION 

Three dimensional seismic data is collected typically 
through the use of vibrators for land data collection and air 
guns for sea data collection. In either event, the data is 
collected, processed, and presented as a three-dimensional volume 
formed of digitized traces. The traces are represented as a 
series of points. each point being described digitally . By 
examining the series of points for significant seismic 
events and by interrelating significant events between traces, 
three dimensional horizons may be developed. 



01 81 21 6 

5 

The processed "three-dimensional seismic data sections 
traditionally have been interpreted by human experts. A tedious 
and time-consuming task that the interpreter must perform is to 
pick seismic events through the volume and to map the time and 

5 amplitude structure of the events. This invention addresses the 
problems of automatically picking and mapping three-dimensional 
seismic event horizons. 

Figure 1A illustrates a typical processed three-dimensional 
seismic data volume 40. The gradations in shade are 

10 generally proportional to changes in amplitude of the seismic 
signal returned as a result of a seismic source such as a 
vibrator or air gun having been activated. 

Figure IB indicates typical three-dimensional horizons 
50a-50e that lie within three-dimensional seismic volume 40 

15 where time is the vertical dimension. These horizons are the end 
result of this invention and may be displayed as shown on paper, 
or may be displayed in the face of a cathode ray tube, or any 
other display apparatus. The horizons may be viewed from any 
direction, of course, and other information about the horizon, 

20 such as the seismic amplitudes on the surface, may also be 
displayed. Horizon, in the strict geophysical sense, means 
"geological interface". As used in this invention, horizon means 
"interpreter's sketch of a surface through seismic data." The 
wavelengths of the source typically imparted into the earth in 

25 seismic exploration are of such a magnitude as to preclude the 
detection of al 1 "geological interfaces". Therefore, the 

horizons that are referred to in this specification are surfaces 
in the three-dimensional seismic volume 40 that are estimated to 
be everywhere locally parallel to the "geologic interfaces". 

30 The automatic production of the horizons 50a-50e is 

accomplished in four steps: 

Step 1: Determine and describe significant points of change 
along each seismic trace. 
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Step 2: Determine pairs of turnings such that the turnings 
within pairs lie on the same three-dimensional seismic event. 
Step 3: Determine seismic event surfaces (horizons). 
Step 4: Produce physical maps of the seismic event surfaces 
5 (horizons). 

Figure 2 illustrates the LISP machine of this invention, a 
LAMBDA machine made by LISP Machines, Inc. The LAMBDA system 
and the use of the LISP programs in the Appendix, are set 

out in detail in the following manuals of LISP Machine, Inc.: 
10 "Introduction to the LISP System" - 1983 

"ZMACS- Introductory Manual" - 1983 
"ZMAIL Manual" - April 1983; 
"Window System Manual" - August 1983; 
"Lisp Machine Manual" - January 1983; 
15 "System 94 Release Notes" - Undated 

in block form, a LISP processor 11 is connected via bus 17 
(a Texas Instruments NuBus™ type bus) to a Motorola 68000 
processor 12. main memory 13, display controller 14 and system 
diagnostic unit 15. Display controller 14 controls video d,sp1ay 
20 20 which is connected to RS232 interface 18. The video display 
20 may be used for displaying the representation of the 
three-dimensional section 40 and for the display of the horizons 
50a-50e. Keyboard 22 and mouse 21 are connected to video display 
20 for interaction therewith. SMDdisk centre 1 1 er 23 control s 
25 disk drive 25 which Is used for storage as is main memory 13. 

4- ^ +-« r<^?32 Dort 29 which in turn is 
Printer 30 is connected to RS23Z porx 

connected through the system diagnostic unit 15 to LISP 
processor 11. main memory 13 and SMDdisk control 1 er 23 . The 
printer 30 is therefore accessible to a storage unit for 
30 receiving a representation of the horizons 50a-50d for printing 
hard copies of the desired images. 

Turning now to Figure 3. Step 1 as carried out by system 10 
will be described. For the purpose of description. each 
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significant point of change along a seismic trace is a turning. 
The description of a turning includes its position in the 
volume, the amplitude of the seismic trace at that position, and 
Its sign (negative or positive). Trace 51 is represented in an 
5 analog manner. Points 52-58 represent turning points in the 
case where turnings are defined by zero-crossings in the first 
derivative. In fact, trace 51 is represented by a series of 
digital numbers. Each number signifies an amplitude. Turning 
points 52, 54, 57 and 58 have positive amplitudes and points 53, 
10 55, and 56 have negative amplitudes. However, the signs of the 
turning are determined by the (N+l)th derivative of the trace and 
therefore the sign of turnings 54, 55 and 56 are positive while 
the sign of turning points 52. 53, 57 and 58 are negative. 
The turning points are determined using zero-crossings in the Nth 
15 derivative. The turnings have the advantage of distinctly 
marking significant changes such as local extrema. saddle points, 
inflections and other high-order changes. This is in contrast to 
prior art methods of looking at the zero crossings of the original 
trace as significant. Exactly which changes are marked can be 
20 controlled by varying the value of N, a parameter of step I. The 
derivative at each sample point, tO, is approximated as follows: 

1. Let t-, tO, and t+ be the times of three consecutive 

sample points. 

2. Let A (t) be the amplitude of the trace at time t. 

25 3, Then the derivative, D, of the trace at tO is computed 

by: 

D (tO) = (A <t+) - A(t-))/2 

This approximation for computing the derivative is derived 
from the parabola P. passing through (t-. A(t-)>. (tO, A(t0)). 
30 and <t+, A(t+)). D (tO) is the slope of P at to. The solution is 
simple because the sample points are regularly spaced at unit 
intervals, allowing the following property of parabolas to be 
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^'Theorem: Let POO-AX Wc Then the der.vatwe of P at a 
value half way between two values, v and w. Is the slope of the 
line through (v, P<v>> end (w, P(w>). That Is = 
5 P'((v+w>/2> - (P(w)-P(v))/(w-v). 

Higher order-der Wat Ives are computed by iteratWe 
applications of the scheme for dcu.-tln. the derWatWe a, 

shown above. . ^ =4 _ i—-,! 

,„ principle, the first derivative defines turnings at local 
10 extreme and saddles of the seismic trace. The first 
extrema a ^tributes of an important turning 

derivativetherefore captures attributes or a k 

point oulte wen. The second der, vat, vedef 
Inflections and straight segments of the trace |nf|ectIons 
the tern, "turning" Is somewhat of a misnomer, but the Inflect. on. 
15 serve just as well or better than the local extrema and 

15 may serve J becomes larger, smoother and 

saddles for a given trace. As n oec 

more complete three-dlm.ns.Pna. horizons are produced. In 
"ect.cT N-7 has the advantage of marKIng changes or mapping 
a are ordlnerlly not noticed by human Inspect.on This may be 
20 seen by an. Inspection of figures SA-6H where figure 6A 
Illustrates traces with N=0. figure 66 Illustrate, figure 
6 C Illustrates N - 2 , and so on with f 1 gure 6H . I »«*^ ]' 

* fKU cteo 1 may be seen in the appendix 
The implementation of this step ^ ls DETERMINE 

beginning on Page 2 thereof under the title STEP 

25 TURNING POINTS". joints in 

By specifying only those points that are turning points 
th. Nth derivative, data compaction is accomplished. Of course 
a variety of other techniques, Including frequency domain f, Iters, 

way be employed to approximate the derivatives. 
30 Figure 4 illustrates traces 60 and SI having the mutual 

nearest neighbor criterion applied. The mutual nearest neighbor 
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criterion was originally devised and used For solving problems In 
computer vision as set out in the Ph.D dissertation "A 
Computational Theory of Spat I o— Tempora 1 Aggregation for Visual 
Analysis of Objects in Dynamic Environments" by Bruce E. 

5 Flinchbaugh, The Ohio State University, June 1980. 

This mutual nearest neighbor criterion has been applied in 
the past to seismic interpretation relating to two dimensional 
seismic sections. This prior art technique involved the 
two-dimensional lineation of wavelets. Wavelets are defined as 

10 that portion of a wiggle trace between amplitude zero-crossings 
when N=0. A 1 i neat ion refers to a two-dimensional horizon in a 
2-D seismic section. However, a technique is applied herein for 
the first time to a three-dimensional volume 40. That is, the 
mutual nearest neighbor criterion is applied in the X direction 

15 and in the Y direction to each and every turning using Nth 
derivative of every trace making up the three-dimensional volume 
40. 

In Figure 4, turning point A is shown connecting to turning 
point B. That is, the nearest neighbor to point A of trace 60 in 

20 trace 61 is point B. Likewise, the nearest neighbor to point D 
of trace 60 is point E of trace 61. The nearest neighbor to 
point F of trace 60 is point G of trace 61. In a similar manner, 
the nearest neighbor to point B of trace 61 is point A of trace 
60. The nearest neighbor of point E of trace 61 is point D of 

25 trace 60. The nearest neighbor to point G of trace 61 is point F 
of trace 60. The nearest neighbor of point C of trace 60 is 
point S of trace 61, but because point B's nearest neighbor is 
not point C, point C does not have a mutual nearest neighbor, 
thus point C is not paired with any turning. It should be noted 

30 that while these traces are shown in a 2-D analog representation, 
the mutual nearest neighbor is based only on the time dimension 
which is the vertical dimension as shown in Fig. 4. 

Turning point 0 has been paired with turning point E. 
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. h- r n». a radically different character from turning 
Turning point E has are d.caly 

P °' nt pen™ I v.. Uke„„e. turning point 0-s 

5 ^Figure 5 nlustrates that turning Points 0 and E mey not be 
Figure Din prohibition against 

Unked nor may turning points F and that the character of 

suc h link, is made based on the assumption that the cha 
. three-dlmenslonal selsm.c ^'^/^^JI^ ™. 

Therefore, the cr.terlon n _„. t ne , gh t,or and homogeneous 

The implementation of - ^ ^ ^ 3 under the 

heading 6 to the heading 

15 HORIZON', extending to page 

■APPLY-U^ATURE^™.^ impl ement at ion. has now 

Each turning point, * hr ^ other turnin gs, 

20 " Ta'pert ^.rrr^ature of the oath formed by the 
, red turnings In successive digital traces Is assessed In both 
Ze l and Y directions. In Figure 7 line segments are drawn 
the X ana Me(Brvt D ath 70 through turnings H, 

between paired turnings to represent path ^ 

^ m The curvature of the patn ax r\ 
25 J. K. L. and «. The curve ^ o ^ K _ ^ 

considering the slope, ml. of the segme 
the slope, mZ. of the segment between K and If 
d ,ffer by »,e than a predetermined value the path 7 ^ 
at K by eliminating -the pe.rs ». « -J^ 
3„ Predetermined value Is a parameter^ * the^ ^ ^ ^ 

"th/Cn,:;; - u - ^;;r^rrr- 

an incident pair in the same direction. Is less 
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curvature. Restricting pairs of turnings in this manner helps to 
rule out implausible three-dimensional seismic horizons. The 
implementation of the low curvature criterion may be studied in 
the Appendix beginning at page 6 under the heading 
5 "APPLY-LOW-CURVATURE-CRITERION". 

Step 2 also includes a requirement of three-dimensional 
continuity. This criterion requires that a pair can only occur 
wherever the pair is a link in a closed loop of turnings. 
Restricting pairs of turnings in this manner helps to assure that 
10 pairs are on seismic events that have three-dimensional extent as 
a horizon. This criterion is based on the assumption that if two 
turnings lie on a three-dimensional seismic event, then other 
turnings lying on the same event surface are nearby. 
Specifically, this criterion is described as follows: 
15 Let Tl and T2 be sets of turnings such that Tl contains the 

turnings from one seismic trace and T2 contains the turnings 
from an adjacent trace In any direction within the seismic data 
section. Let a closed loop of turnings be an ordered list of at 
least three turnings, (to, tl , t2 , ...tn), such that the ti are 
20 distinct and the following pairs exist: 

((t0, tl), (tl, t2) ...(tn-1, tn),(tn, t0)). Thentl in Tl 
and t2 in T2 may only form a pair when a closed loop of turnings 
exists that include tl and t2. 

Closed loops spanning just a few (e.g, four) traces in 
25 closely spaced three-dimensional seismic data can be assumed to 
bound a region of a surface that intersects those traces. 
Reference to Figures 8A and 8B should now be made where, in 
Figure 8A, four traces Ta , Tb, Tc and Td are shown. Turning 
point pairs (a, b) (b, c) (c. d) and (d, a) are shown forming a 
30 closed loop. The four turning points, a-d, form four pairs of 
turnings as indicated and also form a closed loop. On the other 
hand, pairs of turnings (g. e) and (e, f) do not form a closed 
loop and therefore do not define any surface. 
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A top y.ew of a Cosed loop concept is shown ,n figure «• 
, tc anbe seen that area, connected together are areas of four 
turnings f or„„ S four pairs or turnings In a -'•— '«> 
configuration. However, the pairs of turn Ings In O, <m n ,P, 

S p> and <m. P> do not form a dosed loop. « • ^Tj'^ 
pairing between turning points m and p. Therefore, t P 
~ „" is removed and only dlst.net horizons 73- and 74 are 
l l „,ned As may be seen In the Appendix beginmng on page 7 
u^dlr • APP.V-3-0^ N T, N U,TV-CR,T ER .O N -. -n that 1 mp 1 ementet 1 on 
under a^li ^ „,„ of turnings provide 

10 .our turnings forming four pairs of ^ 
three-dimensional continuity. A pair that 

criterion is "isolated" and rejected hQr , zons nave been 

After Steps 1 and 2 are completed, the 
. „ .. h (n the three-dimensional volume 40. Toru.iy 
determined w,thin the tnr rf ^ 

15 identify the horizon, however. entry is 

r^mts defining that horizon and then an 
turning points derini g d1rec tly or indirectly 

remaining turning points that are Reference 

4. ^ +-« the entry turning point are collected, 
connected to the entry , , , ostrate th e criterion as a part 

Should be mad. to Figure 10 a2 
20 of step 3. Turning point d2 Is pa.r 

which in turn ,s Pa.red with C and 
paired w.th turn, ng point c2. Theref . there fore 
d2 are In the same horizon. The op 
er turnings ^n a co ect,on are ion ^ 

25 indirectly) by pairs. Th. tnre e-dlmenslonal 
determination of the informational content and 

eXten Thfs TZZZ ,s described as follows: Let the relation 
, , ? t. tj' indicate that t, and tj are turnings of the same 
(pair) <t., tj). <n transittve . tl an d ta belong to the 

30 hor.zon. This relat on ^ ^ tM 

same horizon when pa r «> and P ^ 

horizon, because pair (t. «> re ,at.on. The specific 

property Is used to complete the pair 
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implementation may be studied in the Appendix beginning on page 7 
under "STEP 3: EXTRACT HORIZONS'. This reformatting results in 
an array format. The array format is appropriate for use in the 
printer 30 and video display 20 of Figure 2 for displaying the 

5 recommended horizons. 

Also, It may be that one horizon is connected to another 
horizon. Step 3 eliminates connections between horizons. A 
horizon can contain only turnings such that no turnings are 
positioned with one directly above the other. This criterion 

10 prevents spiral ing horizons and assists in determining separate 
horizons as follows: 

Let the horizontal position of a turning, ti, be (xi , yi ) , 
Then turnings tl and t2 can only be part of the same horizon if 
(xl, yi) is not (x2, y2), i.e., ff at most one turning of the 

15 horizon lies in any given trace in the volume 40. Thus, each 
turning of a horizon is guaranteed to be "vertically unique". 
Implementation of vertical uniqueness is shown in the Appendix, 
on page 8 under "ADD— TURN I NG-TO-HOR 1 20N ' . 

Figure 9 illustrates the application of this criterion. 

20 That is, a complete loop is not achieved because points al and el 
appear on the same trace Ta. A collection of turnings of a 
horizon may include al and dl, but turning el may not be 
included in a horizon that includes turning points al, bl , cl and 
dl . 

25 The apparatus for both printing and CRT display are well 

known, as are the manner in which that apparatus operates. For a 
specific implementation, please refer to the Appendix beginning 
at page 8 under "STEP 4; DISPLAY HORIZONS". 

While this invention has been described in specific steps 

30 and in specific hardware, it is contemplated that those skilled 
in the art may readily substitute hardware and method steps 
without departing from the scope of this invention which is 
limited only by the appended claims. 
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. .Mode :Lisp; Package : Amos ; Base: 10. r*~ 

Bruce Flinchbaugh 

;AM0S 



* 1===== " == A m7s7 Automatic Mapping of Seismic Events 
Terminology 

a 4« this program, refers to a surface 
The term. Horizon, as used m £^/ documents and im plementations 
computed by the ^^f^.^ Reflector Parallels. Seismic Event 
the surfaces are called, varaously R ^ ^ simplicity; 

Surfaces, and Seismic Events " actual horiz0 ns 

the surfaces do not necessarily correspona 

(geologic interfaces). === ================ === === ====== ' * 

(defflavor AMOS 

'■3D array of seismic amplitodes. indices, Time. and 
S'array.l". M partial derivative of Dana ...... Time. 

nh" order o f the pertiai derivative in Dn-D.ta. 

t„e m*— -IK— I"" 1 «■"-«• iD > " ,tiZO " ' iiP - 
array ovists of Tumin e Point., indices: X and V. 
^cfnorrL. coated b, the Amos method. 

(Horizons nil) are displayed. 

-Window in which maps of the horizons are disp y 

(Window terminal-io) times for maps. 

..List of numbers defining bins of Data times for m p 

^rt'Tnu^ers defining bins of Data amplitudes . 
"•The bins are computed during display (Step 4) . 
(Amplitude -Bins nil) 

) 

•.settable-instance-variables 



(: documentation „ t v n j t 0 use the method 

^The Amos flavor coordinates the Amos «thod. To ^ 

TOke an instance. ^-JjJ^^SoS Message . M , D . and Window 
send the instance the .COMfuic nuiva 
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are optional Init Options." 
)) 



(def method (amos : COMPUTE-HORIZONS) () 

"Constructs the Horizons representation using the Amos method and 
displays maps of the horizons." 
;Step (1). 

(send self rdetermine-turning-points) 
;Step (2) . 

(send self :determine-pairs-of -turnings-on-the-same-horizon) 
;Step (3) . 

(send self : extract-horizons) 
;Step (4). 

(send self rdisplay-horizons)) 



TURNING POINT REPRESENTATION 

(defflavor TURNING 
( 

; Integer, the X coordinate of the location in Dn-Data. 
(X 0) 

; Integer, the Y coordinate. 
(Y 0) 

: Number, the interpolated time of the zero-crossing. 
(Time 0.0) 

; Number, the interpolated amplitude. 
(Amplitude 0.0) 

;Number. the sign (+/-1) of the partial of Dn-Data at Time. 
(Sign 0) 

-.Nearest neighboring turnings in the Time directions ( + /-).. 
(Nt* nil) 
(Nt- nil) 

; Nearest neighboring turnings of the same horizon in the 
;X and Y directions (+/-) . Initially each is just the 
;nearest neighboring turning on the adjacent trace. At the 
;end of Step (2) they connect only paired turnings. 
(Nx+ nil) 
(Nx- nil) 
(Ny+ nil) 
<Ny- nil) 

; Array representing the horizon, if any. of this turning. 
(Horizon nil) 
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) 

0 

-.settable- instance- variables 

(: init-keywords (: specif ication) ) 

(•.documentation . 
"The Turning flavor represents a Turning Point in the Amos method. 
The Specification Init Option. (Time Data Dn-Data) . specifies that 
a turning occurs on the trace between Time and Time+1. Amplitudes 
interpolated for the turnings are used when displaying horizons.")) 

(def method (turning : AFTER :INIT) (init-options) 

"Initializes the Time, Amplitude, and Sign of the turning." 
(let* ((tl (first (get init-options : specif ication))) 

(data (second (get init-options : specif ication)) ) 
(dn-data (third (get init-options -.specification))) 
(t2 (+ tl D) 

(vl (aref dn-data X Y tl)) :Nth partials at tl and t2. 
(v2 (aref dn-data X Y t2) ) 

(al (aref data X Y tl)) ;Amplitudes at tl and t2. 
(a2 (aref data X Y t2))) 
(setq Time (// (- (* (- v2 vl) t2) v2) (- v2 vl)) 

Amplitude (let ((m (// (- a2 al) (- t2 tl)))) 

(+ (* m Time) (- a2 (* m t2)))) 
Sign (if (> v2 vl) 1 -1)))) 



*' STEP 1: DETERMINE TURNING POINTS ^ ^ 

(def method (amos : DETERMINE-TURNING-POINTS) () 
"Constructs the Turnings representation." 
(send self :compute-nth-partial-derivative-of -data) 
(send self :mark-zero-crossings-in-nth-derivative) ) 

(defmethod (amos :COMPUTE-NTH-PARTIAL-DERIVATI VE-OF-DATA) () 

"Computes the nth partial derivative of Data with respect to Time, 
leaving the result in Dn-Data." 

(setq Dn-Data (make-array (array-dimensions Data) :type :art-q» 
:Dn-Data is just Data in case N = 0. 
(copy -array -contents Data Dn-Data) 

(loop, with data* = (make-array (array-dimensions Data) -.type :art-q) 
repeat N do 

(copy-array-contents Dn-Data data*) 

(send self : compute-partial-derivative data* Dn-Data))) 
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(defmethod (amos : COMPUTE-PARTIAL-DERIVATIVE) (data* dl-data*) 
"Computes the partial derivative of DATA* with respect to Tiipe , 
leaving the result in Dl-DATA* . 11 

(loop with (xdim ydim tdim) = (array-dimensions data*) 
for x below xdim do 
(loop for y below ydim do 

; Compute derivatives at times 1 through TDIM - 2: 
(loop for time from 1 to (- tdim 2) do 

(aset (// (- (aref Data x y (+" time 1)) 
(aref Data x y (- time 1))) 
2.0) 

dl-data* x y time)) 
; (Special cases) : 

; Compute derivatives at times 0 and TDIM - 1. 
(aset (- (aref Data x y 1) (aref Data x y 0)) 

dl-data* x y 0) 
(aset (- (aref Data x y (- tdim 1)) 
(aref Data x y (- tdim 2))) 

dl-data* x y (- tdim 1))))) 

(defmethod (amos :MARK-ZERO-CROSSINGS-IN-NTH-DERIVATIVE) () 

"Constructs the Turnings representation, creating one Turning for 
each zero-crossing in Dn-Data. The original implementation of the 
Amos method also created Turnings for zero-touchings ; the touchings 
are omitted here for simplicity." 

(loop with (xdim ydim tdim) - (array-dimensions Dn-Data) 
initially (setq Turnings 

(make-array (list xdim ydim) :type 'art-q)) 
for x below xdim do 
(loop for y below ydim do 

(aset ;Check for turning in each interval [time , time* 1] : 
(loop for time from 0 to (- tdim 2) 

for value-at-time = (aref Dn-Data x y time) 
for value-at-time+1 =- (aref Dn-Data x y 

(+ time 1)) 
if (or (and (plusp value-at-time) 

(minusp value-at-time+1) ) 
(and (minusp value-at-time) 

(plusp value-at-time+1)) ) 
collect (make- instance 

•turning :x x :y y 
.•specification 
■(.time .Data t Dn-Data))) 
Turnings x y)) )) 
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STEP 2: DETERMINE PAIRS OF TURNINGS ON THE SAME HORIZON 

(def method (amos :DETERMINE-PAIRS-OF-TURNINGS-ON-THE-SAME-HORIZON) () 
"Scans the Turnings representation to set the Nx+/- and Nv+/- 
pointers to adjacent turnings of the same horizon. " 
(send self : initialize-nearest-neighbors) 
(send self :apply-pairing-criteria)) 

(def method (amos : INITIALIZE-NEAREST-NEIGHBORS) () 

"Sets the Nt + /- . Nx+/- . and Ny+/- pointers to the nearest 
neighboring turnings in the indicated directions." 
initialize the NW- pointers (these simplify operations later): 
(loop with (xdim ydim) = (array-dimensions Turnings) 
for x below xdim do 
(loop for y below ydim do 

(loop with previous-turning = (first (aref Turnings x y)) 
for current-turning in (cdr (aref Turnings x y)) do 
(send previous-turning :set-nt+ current -turning) 
(send current-turning :set-nt- previous-turning) 
(setq previous-turning current-turning)))) 
initialize the Nx+/- pointers: 

(loop with (xdim ydim) = (array-dimensions Turnings) 
for y below ydim do 
(loop 

for x from 1 below xdim 

for previous-turning-list = (aref Turnings (- x 1) y) 

for current-turning-list = (aref Turnings x y) do 

(send self : initialize-nearest-neighbors-between-two-traces 

:set-nx+ previous-turning-list current-turning-list) 
(send self : initialize-nearest-neighbors-between-two-traces 
:set-nx- current-turning- list previous-turaing-list) )) 
initialize the Ny+/- pointers: 

(loop with (xdim ydim) = (array-dimensions Turnings) 
for x below xdim do 
(loop 

for y from 1 below ydim 

for previous-tuming-list - (aref Turnings x (- y D) 

for current-turning-list = (aref Turnings x y) do 

(send self : initialize-nearest-neighbors-between-two-traces 

:set-ny+ previous-turning-list current -turning- list) 
(send self : initialize-nearest-neighbors-between-two-traces 
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:set-ny- current-turning-list, previous-turning-list))) 

) 

(def method (amos : INITIALIZE-NEAREST-HEIGHBORS-BETWEEN-TWO-TRACES) 
(set-n* turning-list-a turning-list-b) 
-Sets the Nx + /- or Ny+/- pointer, as indicated by SET-N* (e.g. 
:SET-KO. of each turning in the A list to the nearest B list 
turning." 

(loop for a-turning in turning-list-a 

if (not (null turning-list-b)) _ 

do (loop with previous-b-turning = (first turning-list-b) 
for current-b-turning in (cdr turning- list-b) 
until (> (send current-b-turning :time) 

(send a-turning :time)) 
do (setq previous-b-turning current-b-turning) 

finally 

(if (or (null current-b-turning) ;=> 1 b-turnmg. 
(< (abs (- (send a-turning :time) 

(send previous-b-turning 
:time))) 
(abs (- (send a-turning :time) 
(send current-b-turning 
:time))))) 

(send a-tuming set-n» previous-b-turning) 
(send a-turning set-n* current -b- turn ing)))) ) 

(defmethod (amos : APPLY-PAIRING-CRITERIA) () 

"The four criteria in Step 2 of the Amos method are applied to prune 
and adjust the pointers." 

"apply-mutual-nearest-neighbor-and-homogeneous-signs-criteria) 

(send self : apply-low-curvature-criterion) 
(send self :apply-3d-continuity-cri tenon)) 

^(amos^APPLY -MUTUAL-NEAREST-NEIGHBOR- AMD-HOMOGENEOUS-SIGNS-CRITERIA) 

"Scans all turnings, adjusting pointers according to the MNN and HS 
criteria in all four directions." 

Uend ^dlLrmine-nearest-neighbor^ 

(send self : prune- according-to-mutual-nearest-neighbor-criterion) ) 
^farof-DETERMINE-NEAREST-NEIGHBORIHG-CANDIDATES-WITH-EQUAL-SIGNS) () 
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"Scans all turnings. Each turning votes. for the most appropriate 
pairing, according to the MNN and HS criteria, in each direction." 
(loop with (xdim ydim) » (array-dimensions Turnings) 
for x below xdim do 
(loop for y below ydim do 

(loop for turning in (aref Turnings x y) do 

(send self : vote-f or-nearest-equal-sign-candidate 

turning :nx+ :set-nx+) . . 
(send self : vote-f or-nearest-equal-sign-candidate 

turning :nx- :set-nx-) 
(send self : vote-f or-nearest-equal-sign-candidate 

turning :ny+ :set-ny+) 
(send self : vote-f or-nearest-equal- sign-candidate 
turning :ny- :set-ny-))))) 

(def method (amos : VOTE- FOR-NEAREST-EQUAL-SIGN-C AND IDATE) 
(turning n* set-n*) 
"The N* pointer of TURNING is adjusted to the nearest neighboring 
turning in that direction that agrees in sign. An additional 
criterion is applied to preserve the property of mutual nearest 
neighbor pairings that no pairings cross; pointers that may result 
in such a crossing are set to NIL." 
(let ((n* -turning (send turning n*))) 
(if (not (null n* -turning)) 

;If the turnings agree in sign 
;then N* is already the right vote. 

(if (not (- (send turning :sign) (send n* -turning :sign))) 
; Consider changing N* to the next nearest turning 
; without crossing: 

(if (> (send turning :time) (send n*-turning :time)) 
; TURNING is below N* -TURNING so consider the next 
; lowest one: 

(let ((n*t+ (send n* -turning :nt+))) 
(if (and (not (null n*t+)) 

(= (send turning :sign) (send n*t+ :sign))) 
(send turning set-n* n*t+) 
(send turning set-n* nil))) 
; TURNING is above N*-TURNING so consider the next 
;highest one : 

(let ((n*t- (send n*-turning :nt-))) 
(if (and (not (null n*t-)) 

(= (send turning :sign) (send n*t- :sign))) 
(send turning set-n* n*t-) 
(send turning set-n* nil))) 
;N.B. The case of TURNING having the same time as 
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;N*-TURNING is treated as the "above" case.. It would be 
;better to pick the nearest of N*T* and N*T- , as in the 
; original implementation . 
))))) 

(def method 

(aroos : PRUNE- ACCORDING-TO -MUTUAL-NEAREST-NEIGHBOR-CRITERION) () 
"Sets NX/Y/+/- pointers to NIL if they are not in the mutual nearest 
neighbor relation." 

(loop with (xdim ydim) = (array-dimensions Turnings) 
for x from 0 below xdim do 
(loop for y from 0 below ydim do 

(loop for turning in (aref Turnings x y) do 

(send self :prune-according-to-mnn-directionally 

turning :nx+ :set-nx+ :nx-) 
(send self :prune-according-to-mnn-directionally 

turning :nx- :set-nx-' :nx+) 
(send self :prune-according-to-mnn-directionally 

turning :ny+ :set-ny+ :ny-) 
(send self :prune-according-to-mnn-directionally 
turning :ny- :set-ny- :ny+))))) 

(defmethod (amos : PRUNE- ACCORDING -TO-MNN-DIRECTIONALLY) 
(turning n*a set-n*a n*b) 
"Prunes the N*A pointer of TURNING if it does not point to a turning 
in the mutual nearest neighbor relation with TURNING." 
(if (and (not (null (send turning n*a))) 

(not (eq turning (send (send turning n*a) n*b)))) 
(send turning set-n*a nil))) 

(defmethod (amos : APPLY-LOW-CURVATURE-CRITERION) () 

"Scans all pairings, removing those that violate the low curvature 
criterion. " 

(loop with (xdim ydim) = (array-dimensions Turnings) 
for x from 0 below xdim do 
(loop for y from 0 below ydim do 

(loop for turning in (aref Turnings x y) do 
(send self 

; prune -according- to- low-curvature -directionally 
turning :nx- :set-nx- :nx+ :set-nx+) 
(send self 

: prune -according-to- low- curvature -direct ionally 
turning :ny- :set-ny- :ny+ : set-nyO ) ) ) ) 



(defmethod (amos : PRUNE- ACCORDING-TO-LOW-CURVATURE-DIRECTIONALLY) 
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(turning. n»- set-n*- n* + set-n**) 
"Prunes the incident pairings in the * direction if the change in 
dip between the two pairings is greater than D. Dip is computed 
in the (X. Y, Time) space where traces are spaced at unit intervals 
in the X and Y directions and Time is measured in sample points . 
It would be more appropriate to compute dip using the actual trace 
intervals (in meters) and the true times (in milliseconds) . " 
(let ((turning-n*- (send turning n*-)) 
(turning-n*+ (send turning n*+))) 
(cond 

((and (not (null turning-n*-)) 
(not (null turning-n*+)) 
(> (abs (- (- (send turning :time) 

(send turning-n*- :time)) 
(- (send turning-n* * :time) 
(send turning :time)))) 

D)) 

(send turning-n*- set-n*+ nil) 
(send turning set-n*- nil) 
(send turning set-n*+ nil) 

(send turning-n*+ set-n*- nil))))) 

(defmethod (amos : APPLY-3D-C0NTINUITY-CRITERI0N) () 

"Scans all pairings, removing those that violate the three- 
dimensional continuity criterion." 

(loop with (xdim ydim) = (array-dimensions Turnings) 
for x from 0 below xdim do 
(loop for y from 0 below ydim do 

(loop for turning in (aref Turnings x y) do 
(send self 

: prune -according-to-3d-continuity-directionally 
turning :nx+ :set-nx+ :nx- :set-nx- :ny+ :ny-) 
(send self 

: prune-according- to-3d-continuity-directionally 
turning :ny+ :set-ny+ :ny- :set-ny- :nx+ :nx-)) 

))) 

(defmethod (amos : PRUNE- ACC0RDING-T0-3D-C0NTINUITY-DIRECTI0N ALLY) 
(turning na+ set-na+ na- set-na- nb+ nb-) 
"Eliminates the NA+ pairing of TURNING if the pairing is not part of 
a local loop on either side (in the B+ or B- direction) of the 
pairing . " 
(cond 

((and (not (null (send turning na+))) 

(not (or (send self : local-loop-ties-p 
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turning . na+ nb+ .na- nb-) 
(send self : local-loop-ties-p 

turning na* nb- na- nb+)))) 
(send (send turning na+) set-na- nil) 
(send turning set-na+ nil)))) 

(def method (amos : LOCAL-LOOP-TIES-P) (turning &rest n*i) 

"Returns T iff TURNING is reached by following the four n*i 
pointers, starting from TURNING." 
(loop for n* in n*i 

for current-turning first (send turning n*) 

then (send current -turning n*) 
until (null current-turning) 

finally (return (eq turning current -turning) )) ) 



*' ' STEP 3: EXTRACT HORIZONS 



(defmethod (amos : EXTRACT-HORIZONS) 0 m 

"Collects turnings connected by the pairings determined in Step 2 
into horizons, constructing the Horizons representation, 
(loop with (xdim ydim) = (array-dimensions Turnings) 
for x from 0 below xdim do 
(loop for y from 0 below ydim do 

(loop for turning in (aref Turnings x y) 

if (and (null (send turning -.horizon)) 

(or (not (null (send turning :nx+))) 
(not (null (send turning :nx-))) 
(not (null (send turning :ny+))) 
(not (null (send turning :ny-))))) 
do (send self : extract-horizon turning))))) 

(defmethod (amos :EXTRACT-HQRIZON) (turning) + _^ e<3 
"Collects turnings connected to TURNING by the pairings determined 
in Step 2 into a horizon array and adds it to Horizons." 
(let ((horizon (make-array (firstn 2 (array-dimensions Data)) 

:type * art-q))) 
(send self : add- turning-to-horizon turning horizon) 
(push horizon Horizons))) 

(defmethod (amos : ADD-TURNING-TO-HORIZON) (turning horizon) 

"If TURNING is not already part of a horizon and HORIZON does not 
already have a turning at the (X.Y) position of TURNING (enforcing 
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the Vertical Uniqueness criterion),, then .TURNING is added to 
HORIZON and all turnings paired with TURNING are added to HORIZON 
(using the Transitivity of the Pair Relation criterion). 11 
(cond 

((and (not (null turning)) 

(null (send turning : horizon)) 

(null (aref horizon (send turning* :x) (send turning :y)))) 
(send turning : set-horizon horizon) 

(aset turning horizon (send turning :x) (send turning :y)) 
(send self : add- turning-to-horizon (send turning :nx+) horizon) 
(send self : add-turning-to-horizon (send turning :nx-) horizon) 
(send self : add- turning-to-horizon (send turning :ny+) horizon) 
(send self : add-turning-to-horizon (send turning :ny-) horizon)))) 
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« i =««====================== 

STEP 4: DISPLAY HORIZONS 
= _ ================ =========================== S===5S=======!S=====S==,# 

(defmethod (amos : DISPLAY-HORIZONS) (^optional (attribute - time >> 
-Displays color maps of the Horizons, one by one. in Window. The 
ATTRIBUTE values. :TIME or : AMPLITUDE , are indicated by color 
After each display the program waits for any character to be typed 
before proceeding to map the next horizon." 
(if (null Amplitude-Bins) 

(send self : c ompute- amplitude -bins -for-data)) 
(loop for horizon in Horizons do 
(send Window :clear-screen) 

(send self : display-horizon horizon attribute) 
(tyi))) 

;ART-Q-16 of ART-4B-8-1 1 s ; colors 0:15 respectively: 
(defvar 4-BIT-C0L0RS) 

;ART-Q-16 of ART-lB-32-4's ; B&W patterns for colors 0:15: 
(defvar 1-BIT-C0L0RS) 

(defmethod (amos : DISPLAY-HORIZON) (horizon attribute) 

"Displays a color map of the ATTRIBUTE of HORIZON an Window. If 
ATTRIBUTE is : AMPLITUDE all amplitudes lying within the same 
Amplitude-Bin (defined by : COMPUTE- AMPLITUDE-BINS-FOR-DATA) are 
displayed in the same color. If ATTRIBUTE is :TIME the time range 
of the horizon is divided into fourteen equal color bins. The 
f W y=0) corner of the map appears in the lower-left corner of the 
window. The X axis is horizontal and the Y axis is vertical, 
(if (eq attribute :time) 

(send self :compute-time-bins-f or-horizon horizon)) 
(multiple-value-bind (width height) (send window :inside-size) 
(tv:prepare-sheet (window) 

(loop with screen = (tv: sheet-screen-array window) 
and left = (tv: sheet-inside-left window) 
and top «= (tv: sheet-inside-top window) 
and (hxdim hydim) = (array-dimensions horizon) 
with pixel-width = (max 1 (// width hxdim)) 

and pixel-height = (max 1 (// height hydim)) 
with wxmax = (+ left (* pixel-width 

(- (min hxdim width) 1))) 
and wymax = <♦ top (* pixel-height 

(- (min hydim height) 1))) 
with colors = (if (eq 'art- lb (array-type screen)) 
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1-bit-colors 4-bit-colors) 
with alu = tv:alu-ior 

for vry* from top to wymax by pixel-height 
for wy from wymax by (- pixel-height) 
for hy from 0 by 1 do 

(loop for wx from left to wxraax by pixel-width 
for hx from 0 by 1 
for turning = (aref horizon hx hy) 
if (not (null turning)) 

do (bitblt alu pixel-width* pixel-height 
(aref colors 

(send self :color-f or-turning 
turning attribute)) 

0 0 

screen wx wy)))))) 

(def method (amos : COMPUTE- AMPLITUDE-BINS -FOR-DATA) () 

"Sets Amplitude-Bins to a list of numbers defining fourteen 
divisions of the amplitude range of Data. The first seven bans 
represent negative amplitudes and the other seven span the postive 

amplitudes." . 
(loop with (xdim ydim tdim) = (array-dimensions Data) 
and min-amplitude = le34 
and max-amplitude = -ie34 
for x from 0 below xdim do 
(loop for y from 0 below ydim do 

(loop for time from 0 below tdim 

for amplitude = (aref Data x y time) do 
(if (< amplitude min- amplitude) 

(setq min-amplitude amplitude)) 
(if (> amplitude max-amplitude) 

(setq max-amplitude amplitude)))) 

finally 

(setq Amplitude -Bins 

( "o(loop with bin- width = (// C- 0 min-amplitude) 6.0) 
repeat 6 

for bin from (+ min-amplitude bin-width) 

by bin- width 
collect bin) 

!c(loop with bin-width = (// (" max-amplitude 0) 6.0) 
repeat 6 

for bin from (* 0 bin-width) by bin-width 
collect bin) 



Ie34))>) 
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(def method (amos : COMPUTE-TIME-BINS-FOR-HORIZON) (horizon) 
"Sets Time-Bins to a list of numbers defining fourteen equal 
divisions of the time range of HORIZON . Time values are negated so 
that time zero is assigned the largest color number . Thus the time 
color maps are analogous to elevation maps: structural highs of 
horizons are assigned larger color numbers than structural lows." 
(loop with (xdim ydim) = (array-dimensions horizon) 
and min-time = le34 
and max- time = -le34 
for x from 0 below xdim do 
(loop for y from 0 below ydim 

for turning = (aref horizon x y) 
if (not (null turning)) 

do (let ((time (- (send turning :time)))) 

(if (< time min-time) (setq min-time time)) 
(if (> time max-time) (setq max-time time)))) 

finally 

(setq Time-Bins 
•(-le34 

,C(loop with bin-width = (// (- max-time min-time) 13,0) 
repeat 13 

for bin from (+ min-time bin-width) by bin-width 
collect bin) 

le34)))) 

(def method (amos :C0L0R-F0R-TURNING) (turning attribute) 

"Returns a color (bin number) for the ATTRIBUTE of TURNING . " 
(multiple- value-bind 

(tO tl t2 t3 t4 t5 t6 t7 t8 t9 tlO til tl2 tl3 tl4) 
(values-list (selectq attribute 

(: amplitude Amplitude-Bins) 
(:time Time-Bins))) 
(let ((attribute-value (send turning attribute))) 
(if (eq attribute :time) 

; Negate time so structural highs have largest color values: 
(setq attribute-value (- attribute- value) ) ) 
; Binary search for the corresponding attribute -value bin: 
(if (>= attribute -value t7) 

(if (>= attribute -value til) 

(if (>= attribute-value tl3) 

(if (>= attribute-value tl4) 15 14) 
(if (>= attribute-value tl2) 13 12)) 
(if (>= attribute-value t9) 
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(if (>= attribute -value tl0)..ll 10) 

(if (>= attribute -value t8) 9 8))) 
(if (>= attribute-value t3) 

(if (>= attribute -value t5) 

(if (>= attribute -value t6) 7 6) 

(if (>= attribute-value t4) 5 4)) 

(if (>= attribute-value tl) 

(if (>= attribute-value t2) 3 2) 

(if (>= attribute-value tO) 1 0))))))) 

;The remaining variables and functions support the 
; :DISPLAY-HORIZONS method. 



(defvar 1-BIT-C0LORS-SPECIFICATI0N 



• ((#o0000 


#o0000 


#oOOOO 


#oOOOO) ; 


0 1 


-bits (all white background) 


(#ollll 


#olll0 


#ollll 


#ollll) ; 


15 


1-bits (nearly all black). 


(#ol011 


#o!110 


#ollll 


#ollll) ; 


14 




(#ol011 


#olllO 


#o0111 


#ollll) ; 


13 




; (#ol011 


#olllO 


#o0111 


#oll01) ; 


12 


(Chosen to be left out.) 


(#ol010 


#olllO 


#o0111 


#oll01) ; 


11 




(#o0010 


#olllO 


#o0111 


#oll01) ; 


10 




(#o0010 


#olllO 


#o0111 


#o0101) 


9 




(#o0010 


#olllO 


#o0111 


#o0100) 


8 




(#o0000 


#olllO 


#o0111 


#o0100) 


;7 




(#o0000 


#o0110 


#o0111 


#o0100) 


6 




(#o0000 


#o0110 


#o0111 


#oOOOO) 


;5 




(#o0000 


#o0110 


#o0110 


#oOOOO) 


;4 




(#o0000 


#o0010 


#o0110 


#oOOOO) 


;3 




(#00000 


#oOOOO 


#o0110 


#oOOOO) 


;2 




(#o0000 


#oOOOO 


#o0010 


#oOOOO) 


;1 


1-bit (nearly all white). 


(#ollll 


#ollll 


#ollll 


#ollll))) 


;16 


1-bits (all black outline). 



(defun CREATE-4-BIT-C0L0RS () 

(setq 4-bit-colors (make-array '(16) :type 'art-q)) 

(loop for color from 0 to 15 

for array = (make-array '(8 1) :type 'art-4b) do 
(loop for i from 0 to 7 do (aset color array i 0)) 
(aset array 4-bit-colors color))) 

(defun CREATE- 1 -BIT-COLORS () 

(setq 1-bit-colors (make-array '(16) :type *art-q)) 
(loop for color from 0 to 15 

for 1-bit-color in 1-bit-colors-specif ication 
for array = (call # 'tv: make -gray 

() 4 () 4 : spread 1-bit-color) 
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do (aset array 1 -bit-colors.. color))) 

(eval-when (load eval) 
(create-4-bit-colors) 
(create- 1 -bit-colors)) 



Convenience Functions for System Use 

_========================«============= I # 



•.List of all instances of AMOS during a session: 
(defvar * AMOS- INSTANCES* nil) 

(defun APPLY -AMOS-METHOD 

(^optional (data nil) (window termmal-io) (n 1) (d 0.5)) 
"Applies the Amos method to DATA, displaying the resulting time 
structure maps in WINDOW . If DATA is NIL. a 3D array of 
hypothetical data is used obtained via CREATE-SEISMIC-DATA-ARRAY . 
(let ((amos (make -instance 'amos 

:data (if (null data) 

(create-seismic-data-array) 

data) 
: window window 
:n n 
:d d))) 

(push amos *amos- instances*) 
(send amos :conrpute-horizons))) 

(defun CREATE-SEISMIC-DATA-ARRAY 

(^optional (x-dim 10) (y-dim 10) (time-dim 20) 
(trace-generation-function #*sin) 

(dip-generation-function # 'dome-generation-f unction)) 
"Returns an array of hypothetical seismic data suitable for Data m 
the AMOS flavor. This function is not part of the Amos method. 
TRACE-GENERATION -FUNCTION is a function of Time that defines the 
form of every trace . DIP-GENERATION-FUNCTION is a function of X 
and Y that defines the offset of the starting point of each trace 

from zero." m . 

(loop with array = (make-array (list x-dim y-dim time-dim) 

:type 'art-q) 

for x below x-dim do 
(loop for y below y-dim do 

(loop for time below time-dim do 

(aset (funcall trace -gene rat ion- function 
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time 

(f uncall dip-generation-function 
x y))) 

array x y time))) 
finally (return array))) 



(defun DOME- GENERATION - FUNCTI ON (x y) 

"Defines a dome-like structure that peaks at (X=4,Y=6). suitable for 
use with CREATE-SEISMIC-DATA-ARRAY . M 
(» -0.1 (sqrt (+ (* (- x 4) (- x 4)) 

(* (- y 6) (- y 6)))))) 



01 81 218 

j*T 34 

CLAIMS: 

1. A system for automatically producing representations of 
three-dimensional surfaces from processed data in the form of 
digital traces, comprising: 

(a) means connected to receive the digital traces for 
5 identifying events in the digital traces; 

(t>) means connected to receive the identified events for 
determining pairs of events that tentatively lie in the same 
three-dimensional surface represented in a first format; 

(c) means responsive to the determined pairs of events for 
10 reformatting at least one three-dimensional surface in a second 
format ; and 

(storage means connected to receive the second format for 
storing a representation of the three-dimensional surface. 

2. The system of Claim I wherein the means for identifying 
15 comorises: 

means for defining the trace as a series of points in time: 
means for taking the Nth derivative of the trace at each of 
rhe points; 

means for identifying the time of the zero-crossings in the 
20 Nth Derivative; and 

means for establishing the amol itude of the trace at the Nth 
derivative zero-crossing time. 

3. The system of Claim 2 wherein the means for taking the 
Nth derivative comprises means for measuring the slope of the 

25 trace at each point as determined by the immediately preceding 
and following points. 

4. The system of Claim 3 wherein the means for establishing 
the amplitude comorises means for interconnecting the amplitudes 
of the point preceding and following the zero crossing and 

30 marking the amplitude at the crossing of the interconnection 
with the time of the zero-crossing. 
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5. The system of Claim d wherein the means for determining 
oairs of events comprises means for determining pairs of 
turnings further comprising: 

means for comparing the parameters of the turnings in a 
5 first trace with the parameters of the turnings in a second 
trace; 

means for comparing the parameters of the turnings in the 
second trace with the parameters of the turnings in the first 
trace; 

10 means for detecting the differences between the parameters 

of the first turnings in the trace and the parameters of the 
turnings in the second trace, and vice versa; 

means for pairing those turnings having the smallest 
differences, thereby designating the mutual nearest neighbor 
15 pairs of turninas; 

means for repeating the comparison and pairings between 
adjacent traces throughout the three-dimensional seismic volume; 
and 

means for adding mutual nearest neighbor pairing 
20 designations to each turning to form a network of paths o^ 
turn i ngs . 

6. The system of Claim 5 wherein the parameters compared 
comor \ se times of turnings. 

7. The system of Claim 6 further comprising means for 
25 comoaring the signs of the (N+l )th derivatives of the first 

trace at the turning points with the signs of the (N-H)th 
derivatives of the second trace, and for preventing the pairing 
of turnings having unlike signs. 
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8. The system of Claim 7 wherein the means for determining 
pairs of events further comprises: 

means for assessing the curvature of each of the paths of 
turnings of the network in both the X and Y directions; 

means for determining whether the curvature exceeds a 
predetermined value for each of the paths of turnings; and 

means for terminating any of the paths at the turning where 
the curvature exceeds the predetermined value. 

9. The system of Claim 8 further comprising: 

means for determining whether a pair of turnings is a link 
in a closed Iood of turnings and for rejecting any pair of 
turnings that does not constitute a link in a closed loop. 

10. The system of Claim 9 wherein the means for assessing 
comprises means for successively measuring the slope (ml) of the 
segment between a turning and a preceding turning, and the slope 
(m2) of the seoment between the turning and a succeeding 
turning. 

11. The system of Claim 10 wherein the means for 
determining comprises means for measuring the difference between 
slopes ml and m2. 

12. The system of Claim 11 wherein the means for 
terminating comprises means for comparing the difference with a 
predetermined value and for rejecting the segments if the 
difference exceeds the predetermined value. 

13. The system of Claim 12 wherein the means for 
reformatting comprises: 

means for accessing any one of the turnings in the network: 
means for accessing every turning linked direct lv or 

indirectly to the first accessed turning; and 

means for establishing at least one parameter for each or 

the accessed turnings, defining a second format. 
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14. The system of Claim 13 further comprising means for 
rejecting any turning during the accessing that has been 
previously accessed. 

15. i he system of Claim 14 wherein at least one parameter 
5 comprises the address of the corresponding turning in the first 

format . 

16. The system of Claim 14 wherein at least one parameter 
comorises the time of the corresponding turning in the first 
format . 

10 17. The system of Claim 16 wherein another parameter 

comprises the amolitude of the corresponding turning in the 
first format. 

L8. The system of Claim 17 further comprising indicia 
means, connected to receive and to manifest a representation of 
15 the second format. 

19. A method of automatically producing representations of 
three-dimensional recommended horizons from processed seismic 
data in the form of digital traces, comprising the steps of: 

(a) identifying turnings in the digital traces; 
20 (b) determining pairs of identified turnings that 

tentatively lie in the same horizon represented in a first 
format ; 

(c) reformatting at least one horizon, represented in a 
second format ; and 
25 (d) storing the second format as a representation of the 

recommended hor i zon . 
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20. The method of Claim 19 wherein the step of identifying 
turnings comorises the steps of: 

defining the trace as a series of equally spaced turnings in 

time; 

taking the Nth derivative of the trace at all of the equally 

spaced turnings; 

identifying the time of the zero-crossings in the Nth 

derivative; 

estaolisning the ampl itudes of the traces at the times of 
tne zero-crossings: and 

storing the time and amplitude of the trace at the 
zero-crossing of the Nth derivative. 



21. The method of Claim 20 wherein the step of determining 
pairs comorises the steps of: 

comparing the parameters of the turnings in a first trace 
with the parameters of the turnings in a second trace; 

comparing the parameters of the turnings in a first trace 
with the parameters of the turnings in a second trace; 

comoaring the parameters of the turnings in the second trace 
with the parameters of the turnings in the first trace; 

detecting the differences between the parameters of the 
turnings of the first trace and the parameters of the turnings 
of tne second trace, and vice versa; 

pairing those turnings having the smallest differences 
thereby designating the mutual nearest neighbor pairs of 
turnings; and 

repeating the above five steps for all of the traces in the 
three-dimensional seismic section to form a network of paths of 
turn i ngs . 
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22. The method of Claim 21 wherein the step of determining 
pairs further comprises the steps of: 

assessing the curvature of each of the paths of turnings of 
the network in both the X and Y direction; 
5 determining whether the curvature exceeds a predetermined 

value for eoch o'f the paths of turnings; and 

terminating any ' of the paths at the turning where the 
curvature exceeds the predetermined value. 

23. The system of Claim 22 wherein the step of reformatting 

10 comprises the steps of: 

accessing any one of the turnings in the network; 

accessing every turning linked directly or indirectly to the 

first accessed turning; and 

establishing at least one parameter for each of the accessed 

15 turnings, defining a second format. 
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